Probing the time dependence of dark energy 
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A new method to investigate a possible time-dependence of the dark energy equation of state w 
is proposed. We apply this methodology to two of the most recent data sets of type la supernova 
(Union2 and SDSS) and the baryon acoustic oscillation peak at z = 0.35. For some combinations of 
these data, we show that there is a clear departure from the standard ACDM model at intermediary 
redshifts, although a non-evolving dark energy component {dw/dz = 0) cannot be ruled out by 
these data. The approach developed here may be useful to probe a possible evolving dark energy 
component when applied to upcoming observational data. 



PACS numbers: 98.80.-k, 95.36.- 



I. INTRODUCTION 
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Approaches to modeling the dark energy are based ei- 
ther on a particular choice of its equation of state (EoS) 
parameter w{z) or on modifications of gravity at very 
large scales. In the context of general rclativistic models, 
at least three different ways may be followed in order to 
find the dark energy EoS from observations. The first 
and most direct one is to solve the scalar field equations 
for a particular theory, which clearly cannot provide a 
model-independent parameter space to be compared with 
the observational data. Other possibilities are either to 
build a functional form for w{z) in terms of its current 
value wq and of its time-dependence w' or to perform 
a parameter- free approach, such as binned EoS, decom- 
position into orthogonal basis and principal component 
analysis. 

However, irrespective of the choice made, it seems clear 
nowadays that it is extremely difficult, if not impossible, 
to single out the best-fit dark energy model on observa- 
tional data basis only. As well observed in Ref. one 
may, instead of trying to find which dark energy scenario 
is correct, test which models can be ruled out by the 
available observations. This, of course, will not unveil 
the nature of the dark energy but may, with the current 
and upcoming observational data, diminish considerably 
the range of possibilities. 

An interesting example in this direction involves two 
of the favorite candidates for dark energy, namely, the 
vacuum energy density (A) and a dynamical scalar field 
($). Among other things, what obscrvationally may dif- 
fer these two candidates for dark energy is that, in the 
former case, the EoS associated with A is constant along 
the evolution of the Universe [w = —1) whereas in generic 
quintessence scenarios w{z) is a function of time 0. 
Therefore, taking this important difference into account, 
one may conclude that if any observable deviation from a 



constant EoS is consistently found, this naturally poses a 
problem for any model based on this assumption, which 
includes our standard cosmological scenario. 

In this paper, we propose a method to investigate a 
possible time-dependence of the dark energy equation of 
state from current observational data. Differently from 
other approaches that use parameterizations of w{z) 0- 
Q or pde(-z) 0j we use only the Taylor expansion of 
the dark energy density pde(-z) around different values 
of z and the conservation equation as a recurrence for- 
mula, so that the pueIz) derivatives can be directly re- 
lated to w{z) and its derivatives. In doing this, we reduce 
considerably the smearing effect (due to the multiple in- 
tegrals that relate 'w{z) to cosmological distances) that 
makes constraining w{z) extremely difficult (see, e.g., 
0). We note that, unhkc the first order parameteri- 
zations of u'(z). in which the time-dependence is deter- 
mined uniquely by the value of w' at a specific z (usually 
z = 0), in this approach such variation may be verified 
at different values of z by changing the expansion center 
z* at small redshift intervals. From this, we build w{z) 
and w'(z) from the data and discuss the behavior of the 
w — z and w' — z spaces. 

We use in our analyses two of the most recent SNe 
la data sets, namely, the updated compilation of Union 
sample (U2) [| and the nearby -h SDSS ESSENCE 
SNLS + Hubble Space Telescope (HST) set of 288 SNe 
la discussed in Ref. m (thr oughout this paper we refer 
to this set as SDSS compilation). We consider two sub- 
samples of this latter compilation that use SALT2 [11 1 
and MLCS2k2 ^ SN la light-curve fitting methods. 
Along with the SNe la data, and to diminish the de- 
generacy between the two free parameters of the model 
we also use the baryonic acoustic oscillation (BAO) peak 
at ZBAO = 0.35 of Ref. 111. 



II. THE METHOD 
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Let us consider a Friedmann-Robertson- Walker uni- 
verse dominated by non-relativistic matter (baryonic and 
dark) and some form of dark energy. The Fricdmann 
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FIG. 1: Contours of Ax^ ~ 6.17 in the plane w{z* = 0) — w'{z* — 0). Solid contours stand for our second order expansion of 
Pde{z) while dotted and dashed ones correspond, respectively, to CPL and the parameterization of Ref. 0]. Black dots indicate 
the ACDM model. 



equation for this scenario is given by 

2 SttG r 



Puia) + PDE(a) 



(1) 



where H and k stand for the Hubble and curvature pa- 
rameters, a — 1/(1 + z) is the cosmological scale factor, 
and pm and pde are, respectively, the non-relativistic 
matter and dark energy density parameters. 

Now, let us suppose that Pbe{z) is an analytic function 
in the range {z* — e_, z* + e+), so that its Taylor series. 



PBEiz) = P-de{z*) + p'de {z ~ z*) + 



(2) 
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converges to puEiz) in this range. Since the components 
of the cosmic fluid are separately conserved, the dark 
energy density pue{z) and its EoS parameter w{z) are 
directly related by the conservation equation 



PBEiz) = 3pde(^) 



1 + w{z) 

1 + 2 



(3) 



Note that we can use the above equation as a recurrence 
formula to write the derivatives of Pde at z = z* in terms 
of the EoS parameters derivatives, i.e.. 
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P'iE 



1Pbe{z*) , 

[j^ + {3 w'{z*) + z*)]pue{z*) , 



where 7 = 3 [1 + w{z*)]/{l + z*) and a prime denotes 
differentiation with respect to z. From the above results, 
we can rewrite the Friedmann equation ([T|) as 



= Hi 



nM,o{i + zf + —ip{z) + nk,Q{i + zf , (4) 



where E = 1 - Om,o - ^k.o, 



ip{z) = l+7(z-z*) + - 
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and the subscript "0" denotes present-day quantities. In 
what follows, we will restrict our analysis up to the sec- 
ond order expansion of pde(2). The practical reason 
for this choice is that it involves the minimal number of 
parameters [w^z) and w'[z)] required to probe the EoS 
time-dependence. However, in doing so it is important to 
observe that we are also assuming that the real function 
describing the dark energy density is a well-behaved func- 
tion of z with a smooth behavior in the z* ± e± interval, 
so that the higher order terms are negligible compared 
to those considered. This assumption in turn is well jus- 
tified since the result w ~ const, seems to be in good 
agreement with most of the current data (see, e.g., [l3|). 



III. APPLICATION, RESULTS AND 
DISCUSSION 



In order to discuss the current observational con- 
straints on w{z) and w'{z) we use two of the most re- 
cent SNe la compilations available, namely, the Union2 
sample of Ref. Q and the two compilations of the SDSS 
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FIG. 2: Contours of Ax^ ~ 6.17 in the plane w — w' for different values of z* . Here, solid, dashed and dotted curves stand for 
the Ax^ contours from U2 + BAG, SDSS (SALT2) + BAG and SDSS (MLCS2k2) + BAG, respectively. 



collaborations discussed in Rcf. [10[. The Union2 sam- 
ple^ is an update of the original Union compilation. It 
comprises 557 data points including recent large sam- 
ples from other surveys and uses SALT2 for SN la light- 
curve fitting. The SDSS compilation^ of 288 SNe la uses 
both SALT2 and MLCS2k2 hght-curve fitters^ and is dis- 
tributed in redshift interval 0.02 < z < 1.55. 

Along with the SNe la data, and to diminish the de- 
generacy between the dark energy parameters f2M,0j w{z) 
and we also use the BAO parameter 

A = Dy^^^^ = 0.469 ± 0.017 , (5) 

^BAO 

where I?y = [r^{zBAo)zBAo/H{zBAo)]^^^ is the so-called 
dilation scale, defined in terms of the comoving distance 
r to Zbao- The present value of the Hubble parameter 
Hq is taken as a nuisance parameter and we marginalize 
over it. In our statistical analysis, therefore, we minimize 



^ http://www.supcrnova.lbl.gov/ 
^ http://www.sdss.org/ 

^ It is worth emphasizing that SALT2 fitter does not provide a 
cosmology-independent distance estimate for each SN, since some 
parameters in the calibration process are determined in a si- 
multaneous fit with cosmological parameters to the Hubble dia- 
gram (see lis! for a discussion). See also [Tol for a discussion on 
possible systematic biases from the multicolor light curve shape 
method. 



the function ~ Asnc + Xbao ' "which takes into account 
all the data sets mentioned above. Motivated by the re- 
cent results of the CMB power spectrum [l^ , we assume 
spatial flatness in the following analyses. 

Figure 1 shows contours of Ax^ = 6.17 arising from 
the joint analysis U2 + BAO (Panel la) , SDSS (SALT2) 
+ BAO (Panel lb) and SDSS (MLCS2k2)+BAO (Panel 
Ic) for the second order expansion around z* = (solid 
contour). For the sake of comparison we also display the 
corresponding Ax^ contours for two different first order 
parameterizations, i.e., 

2; 

W{z) = Wq + Wa— , (6) 

1 + Z 

the so-called CPL parameterization Q (dotted contour) 
and 

W{Z) Wo + lUa——^ 2~ ' 

1 + 2;^ 

proposed in Ref. Q (dashed contour)"'. We note that, 
although there is space for both negative and positive 



Differently from CPL expression, this latter EoS parameteriza- 
tion is a bounded function of the redshift throughout the entire 
cosmic evolution, which allows to study the effects of a time vary- 
ing EoS component to the distant future at z ~ — 1 as well as 
back to the last scattering surface of the CMB at z ~ 1100 (see 
also for other first order parameterizations). 
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FIG. 3: Evolution of w and w' arising from tlie combinations of data discussed in the text. Error bars correspond to Ax^ = 1 
(Icr). Clearly, no evidence for an evolving EoS is found over the entire redshift interval. Note, however, that for SDSS 
(MLCS2k2) + BAO data a clear departure from the ACDM case {w = —1 and w' — 0) is detected for some values of z* . 



values of w'{z* = 0) (or, equivalently, Wa), as well as 
for phantom-like and quintessence-like Q dark en- 
ergy, the allowed parameter space for our second order 
expansion can be considerably reduced relative to the 
one produced by the first order EoS parameterizations^. 
Wc believe that this constraining power shown in these 
analyses may be a consequence of the reduction of the 
smearing effect mentioned earlier Q, since in our ap- 
proach the parameters w(z*) and w'{z*) appear in the 
expansion rate (Eq. ^) not as an exponent of the pde 
function (as in the case of first order parameterizations) . 

However, as mentioned earlier, the most interesting as- 
pect in the above approach is that a possible time evo- 
lution of w{z) and of its derivatives can be verified from 
the data at different values of z* by changing the expan- 
sion center at small redshift intervals. This procedure 
differs fundamentally from usual first order parameteri- 
zations in that variations of the EoS parameter at z ^ 
cannot be detected if w{z) is a smooth function around 
z ~ 0. In other words, this amounts to saying that in 



For a subset of model parameters, both parameterizations can 
be approximated by their second order expansion. As an exam- 
ple, by taking the best-fit values of uiq and uiq from U2 sample 
(see Table I), we found that the relative error Ap/p < 10% at 
2 ~ 1.5, where Ap is the diff'erence between the energy density 
derived from the above parameterizations and their second order 
expansions. 



the approach discussed in this paper, finding the pair 
{w{z* = 0),w'{z* = 0)) = (-1,0) inside a significant 
confidence region cannot be taken as a definitive argu- 
ment for a non-evolving dark energy since a single value 
of w' ^ at any z* ^ is enough to ensure its time- 
dependence. 

This possibility is tested in Figures 2a-2i, where we 
show contours of Ax^ = 6.17 in the w - w' plane for the 
three combinations of data discussed above. The con- 
tours are obtained as follows: 

1. fix the expansion center z* in Eq. (j?]); 

2. perform the statistical analysis described above to 
obtain the values of Om.Oi w and w' and their cor- 
responding error bars at z*; 

3. set z* — !• z* + Az* and repeat the previous step to 
estimate the above parameters at z* -I- Az*. 

In our analysis, we divided the entire interval in Az* = 

0. 2 from z = up to the highest redshift of our sample, 

1. e., z = 1.4. In Panels 2a-2i, solid, dashed and dot- 
ted curves stand for the Ax^ contours from U2 + BAO, 
SDSS (SALT2) + BAO and SDSS (MLCS2k2) + BAO, 
respectively. For all combinations of SNe la + BAO, it is 
clear that the second order expansion is fully compatible 
with a non-evolving EoS {w' = 0) for the entire redshift 
interval considered. The same, however, cannot be said 
when we consider specifically the standard ACDM case 
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TABLE I: Results for w and w' for the combinations of data discussed in the text. The error bars correspond to la (Ax^ ~ 1) 



and the combination SDSS (MLCS2k2) + BAO (dotted 
line). In this case, pde = const, or, equivalently, w = —1 
and If' = (sec Eq. [2]), which is off from the best- fit 
value by > 2(7 (c.l.) for several values of z*. For the sake 
of completeness, we also show the resulting contours for 
two analyses with z* < 0, i.e., z* ~ —0.4 (Panel 2a) and 
z* = -0.2 (Panel 2b). 

Figure 3 shows the quantities w and w' as functions of 
the redshift for the combination of data discussed above. 
The error bars correspond to the confidence intervals ob- 
tained for Ax^ — 1 • Following the curve indicated by the 
best-fit points, the w — z plane for all combinations seems 
to indicate a freezing behavior, which in agreement with 
recent discussions using reconstructing techniques such 
as the principal component analysis jlSj and the maxi- 
mum entropy method [l^. From these Panels, we can 
see more clearly that, although compatible with a time- 
independent EoS parameter (w' = 0), some combinations 
of the current data (SDSS (MLCS2k2) -t- BAO) can show 
a departure from the ACDM case (w = — 1 and w' ~ 0) 
for some values of z*. More quantitatively, this can be 
seen in Table I, where we display the numerical values ob- 
tained in our analyses. For completeness, we also tested 
the influence of a nonvanishing curvature parameter on 
these results (see, e.g., [2^ for a discussion). By con- 
sidering the current interval from WMAP 4- Hq + BAO 
data, i.e., flk = -0.0023l^;^^^6 EH' ^e verified that aU 
the above conclusions remain unchanged. 



IV. CONCLUSIONS 

We have proposed a new method to constrain a pos- 
sible time-dependence of the dark energy from current 



observational data. By expanding the function pde(-z) 
around different values of z and using the conservation 
equation as a recurrence formula, we have shown a clear 
compatibility of current SNe and BAO data with a non- 
evolving EoS for all samples considered. Our method is 
based on the assumption that the dark energy density 
has a smooth behavior up to z ^ 1.5, so that higher 
order terms in the poE expansion are negligible when 
compared with the first and second ones. Besides involv- 
ing the minimal number of parameters w^z) and w'{z) 
required to probe the EoS time-dependence (see Eq. [¥]), 
this seems to be a plausible assumption since a dark en- 
ergy behavior somewhat close to a cosmological constant 
(pde = const.) is fully supported by many analyses (see, 
e-g-, [13 and refs. therein). We have also discussed some 
advantages of this approach over first order parameteri- 
zations commonly used in the literature and shown that 
the so-called smearing effect in the determinations of the 
w{z) evolution may be greatly reduced. We believe that 
the next generation of experiments dedicated to mea- 
sure the dark energy properties (mainly those measur- 
ing the expansion history from high-z SNe la, baryon 
oscillations, etc. - see, e.g., [IJ]) will probe cosmology 
with sufficient accuracy to decide between an evolving 
and a non-evolving dark energy EoS. In this regard, the 
methodology discussed here can be useful. 
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